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Abstract 

We point out that in the context of quantum fields in time dependent classical backgrounds, 
the number of created particles with a given momentum largely deviates about its mean value. 
Since the corresponding Fourier modes are nonlocal, this deviation shows that the expectation 
value of the number operator can only make sense in an ensemble of spacetimes. Using a complete 
orthonormal family of localized wave packets, we show how an ensemble interpretation can be 
given to cosmological particle creation in local terms. The reheating process following inflation is 
reexamined in the light of this construction. 
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One of the main assumptions of quantum theory is that a physical quantity is well defined 
in a state if the state is an eigenvector of the operator corresponding to that quantity. When 
this is not the case, the quantity does not acquire a definite value and one can only talk 
about the probability of seeing an outcome of a measurement which can be understood in 
an ensemble picture. The expectation value of the corresponding operator gives the average 
of the measurements done on the identical copies of the same system in the ensemble. 
The aim of this work is to reinterpret cosmological particle creation in time dependent 
classical backgrounds by taking into account this basic ensemble picture and discuss possible 
implications of this reinterpretation for inflation. 

To illustrate our point, let us start by reviewing the well known quantum mechanical 
harmonic oscillator with externally fixed time dependent frequency 0. This simple system 
essentially mimics the basic features of cosmological particle creation. Similar to the usual 
harmonic oscillator, by introducing suitable time varying ladder operators a and the 
Hamiltonian can be brought into the form H = u a^a + 1/2 (see e.g. P]). Therefore, the 
instantaneous ground state of the system at time t can be defined by a|0( >= and a com- 
plete orthonormal basis of energy eigenvectors can be introduced as \n >= (o-^)"/ Vn^- \0t >. 

Time dependence of the ladder operators can be solved by a Bogoligov transformation 

a = aao + (3* aj, 

a' = f3 ao + a* Oq, (1) 

where Oq and Oq are the constant operators at time to and a, (3 are complex functions of 
time with — |/3p = l. It immediately follows from this construction that the ground 
state of the system at time Iq, i.e. ao|Oo >= 0, will not remain to be the ground state at a 
later time. Instead, the expectation value of the number operator = a^a can be found as 
< N >=< Oo|iV|Oo >= which shows that on the average this state contains quanta 
at time t. 

Since |Oo > is not an eigenvector of the number operator A^, the above average can only 
be understood in an ensemble picture. Therefore, it is crucial to determine the deviations 
about the mean value. If there is a sharp maximum around the mean with a small deviation, 
then to a very good approximation the ensemble interpretation can be ignored. Expanding 
|Oo > in the basis vectors \n > and using ([1]), one can determine the probability P2n of 
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finding 2n number of quanta in |Oo > as 

{2n)\ 



2n 



22«(n!)2 |a|2n+i- 

Unfortunately, this distribution does not have a maximum around the average On the 
contrary, it is a decreasing function of n. This shows that although on the average |Oo > 
contains quanta, the most probable outcome of a single measurement is the ground state 
|0t > with no quanta. Not surprisingly, the deviation AA^ about the average is large. Using 
(AA^)2 =< > _ < j\r >2 Qj^g calculate the relative deviation as 

V2[-^>V2, (3) 



which shows that the statistical fluctuations about the mean is quite big. 

Let us now consider a free, massive, real scalar field on a cosmological Robertson- Walker 
background ds"^ = —dt^ + a'^{dx'^ + dy"^ + dz"^). One may assume that in addition to the scale 
factor a, the mass parameter M may also depend on t so that particle creation can arise 
due to a coupling to an external time varying field as in the case of reheating. Expanding 
the field in Fourier modes, one can see that this free field theory is nothing but an infinite 
collection of harmonic oscillators each of which is labeled by a comoving wave vector k and a 
time dependent frequency ul = + k"^ /a^ — 9H'^/A — 3H/2, where H = d/a is the Hubble 
parameter (see e.g. 2|]). As in the harmonic oscillator, one can introduce time dependent 
ladder operators and obeying [ok, aj^./] = 5(k — k'), and construct the usual Fock space 
on the instantaneous ground state defined by ak\Ot >= 0. 

The particle creation effects in this setup can also be studied by introducing a Bogoligov 
transformation 

Ok = Oikakito) + PlaLi^ito), 

alk = /?fcak(to) + afcaLk(^o), (4) 

where the time dependent complex coefficients satisfy |/5fcP = 1 and the modes having 

momenta ±k are coupled because of the reality of the field. As before the ground state of 
the system at time to, defined by ak(to)|Oo >= 0, becomes a multiparticle state at a later 
time. 

To find the mean number of created particles in the mode k, one can calculate the 
expectation value of the number operator A^k = cit.O'k, which reads 

< A^k >=< Oo|A^k|Oo >= WS{0). (5) 



The singularity in the delta function can be explained as S{0) = V/{2'k)^, where V is the 
total comoving volume of the space. Defining the comoving number density as rik = A^k/^ 
one obtains 

\0 P 

< nk >=< OolriklOo >= (6) 

which is the familiar expression used in the literature (see e.g. [Sj). 

How one should interpret ([6]) in quantum theory? Since the mode k is globally defined 
and |Oo > is not an eigenvector of nk, ([6]) only gives the average of observations carried out in 
an ensemble of spacetimes. Due to the large deviation about the mean (as in the harmonic 
oscillator discussed above), it is meaningless to use the expectation value (jH]) in a single 
realization. The situation is different for the total number density n = J d?kn\^ or for the 
number density operator integrated out in a finite momentum domain as one can see that 
in both cases the relative deviation in the state |Oo > drops like which vanishes in 

the infinite volume limit. However, in the strict limit the momentum modes are completely 
dislocalized and the density concept becomes ill defined. 

It is clear that the above mentioned problem with ([6]) arises due to the nonlocal nature 
of Fourier modes. To reinterpret particle creation in local terms, we introduce a new set 
of ladder operators corresponding to localized wave packets as follows. Let e > be 
an arbitrary comoving momentum scale and introduce two vectors J = (ji,j2;j3) and n = 
(ni,n2,n3) with mie^'er entries. Let 

"J" = ^ exp(-27rm.k/ e) ak, (7) 
where we use a shorthand notation for a three dimensional momentum integral 



'Jie •'32<^ -'Jse 



/ dkidk2dk^. (8) 

J ? .■!£ 



These operators obey [ajn, akm] = ^JK'^nm and it is easy to see that ajn is the annihilation 
operator for the mode function 

/jn(x) = ^^^^^3^^ ^exp(-ik.(x-27rn/e)), (9) 

which is peaked around the comoving position 27rn/e with a spread 1/e in each direction. 
These localized wave packets form a complete orthonormal family and thus a legitimate basis 
in the Hilbert space, as good as the Fourier modes (2-7?)"^/^ exp(2k.x). Roughly speaking, 



4 



(Q can be viewed as a mode function of a particle localized in a cubic volume of side length 
27r/e, carrying an average comoving momentum Je. 

The cosmological production of these excitations can be determined from the expectation 
value of the number operator A'^jn = aj^aj^- The number density rij^ is defined by dividing 
A^jn to the volume of the cube (27r/e)^. The expectation value of njn in the state |Oo > can 
be calculated as 



< njn >=< Oo|njn|Oo >= \(3k\ . (10) 

The e dependence of this expression is hidden in the definition of the momentum integral 
(IE]) and due to the underlying translational invariance it does not depend on n. Contrary 
to (E]), it is possible to interpret ( ITOi) even in a single spacetime as it gives the average of 
independent observations carried out in different spatial cubic regions labeled by n. Namely, 
there is no difference in taking the average over the same cube in different spacetimes or over 
different cubes of the same spacetime. In this picture the locality of the particle creation 
process is evident. 

Consider now the the expectation value of the total number density rin = Z^j'^Jn, which 
determines the average number of excitations in each cubic region. Not surprisingly, < > 
is equal to < n > defined in the paragraph below ([6]), which shows that both pictures agree 
on the average number of created particles. However, instead of decreasing like 1 / VV, the 



relative deviation of rin in the state |Oo > behaves differently as (Ann/ < ""ti >) ~ e'^^^ [5|.^ 
Therefore, even in the infinite volume limit the number of excitations varies from one cubic 
region to another. Although the background is homogeneous, there appears inhomogeneities 
due to the quantum mechanical fluctuations and locality of the particle creation process. 

For small e, the volume of the cubic regions that are observed becomes large and the 
deviation in the number density naturally diminishes. For large e the opposite is true. In 
an experimental setup e should be determined by the details of the measurement. On the 
other hand in an expanding universe the Hubble radius sets a natural scale since by causality 
particle creation effects should take place independently in each Hubble volume and local 
interactions cannot restore homogeneity. 



^ One simple way to see this dependence is to note that in a homogeneous background the number of 
particles in a volume is proportional to and the relative deviation should decrease with the 

square root of the number of degrees of freedom. 
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As an illustrative application of the above formalism, let us finally consider the reheating 
process after inflation in the parametric resonance regime. Assume that particles are pro- 
duced around a momentum A; in a band having a width Ak with equal average production. 
For each wave packet in the band, the deviation in the number of created quanta is approx- 
imately given by the quantum mechanical result (jSj). Since the momenta are split up by 
e, there are nearly k'^Ak/e^ number of different wave packets in this band. Therefore, the 
relative deviation of the total number of created particles should behave like (fc^Afc/e^) 



(note the e dependence). As shown in [^l, in a simple chaotic model with the infiaton po- 
tential m^0^, scalar x particles can be produced in the parametric resonance regime due to 



the coupling —gcpx- Iii this model, Ak ^ k ^ gm(f)o, where 0o is the initial infiaton 



amplitude, and t 
infiaton mass m 



le Hubble radius at the end of the resonance is approximately given by the 
6j. Setting e ~ m, the order of the relative deviation in the number density 
in each Hubble volume becomes roughly {g(j)Q/m)~^/^. It turns out that a rigorous deriva- 
tion, which takes into account the details like the cosmological expansion, agrees with this 
naive estimate j^. For a realistic set of values given in j^, namely for g = 10^^, (pQ = Mp/2 
and m = lO^^Mp, the order of the relative deviation becomes 10~^. This indicates that 
the "perfect" homogeneity attained after inflation, which is better than one part in e^°, is 
brokendown on Hubble length scales during reheating. 

Acknowledgments 

This work is partially supported by Turkish Academy of Sciences via Young Investigator 
Award Program (TUBA-GESiP). 



[1] R. M. Wald, Quantum field theory in curved space-time and black hole thermodynamics, 

Chicago, USA: Univ. Pr. (1994). 
[2] S. A. Fulling, Aspects of quantum field theory in curved space-time, London Math. Soc. Student 

Texts 17 (1989) 1. 

[3] N. D. Birrell and P. C. W. Davies, Quantum fields in curved space, Cambridge, Uk: Univ. Pr. 
(1982). 

[4] S. W. Hawking, Particle creation by black holes, Commun. Math. Phys. 43 (1975) 199. 

6 



[5] A. Kaya, Quantum mechanical breakdown of perfect homogeneity in reheating after inflation, 

larXiv:0802.2T88] [hep-th]. 
[6] L. Kofman, A. D. Linde and A. A. Starobinsky, Towards the theory of reheating after inflation, 

Phys. Rev. D 56 (1997) 3258, ihep-ph/9704452 



7 



